New structure around 3250 MeV in the baryonic B decay and the D*(2400)iV molecular hadron 
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In this work, we first propose the isovector «D*,(2400)° molecular state to explain the enhancement structure 
around 3250 MeV (X r (3250)°) in the E ++ 7r~7r~ invariant mass spectrum newly observed by the BaBar Collab- 
oration. Under this molecular state configuration, both the analysis of the mass spectrum and the study of its 
dominant decay channel can well depict its resonance parameters measured by BaBar. Our investigation also 
shows that the isovector ;i£)q(2400)° molecular state can decay into I.+ + n~n~, which is consistent with the ex- 
perimental observation. These studies provide the direct support to the isovector nD* Q (24QQ)° molecular state 
assignment to X c (3250)°. 

PACS numbers: 14.20.Pt, 14.20.Lq 



The BaBar Collaboration recently reported a new enhance- 
ment structure in the X ++ 7r~7r~ invariant mass spectrum of the 
B — > Y,^pn~n~ decay. The mass and width are given by 
M = 3245 + 20 MeV and T = 108 ± 6 MeV, respectively 
[1]. We will call the new structure as X c (3250)° in this work. 
In terms of the decay channel observed, we conclude that 
X ( .(3250)° is an isotriplet with charm number C — +1. It is an 
intriguing research topic to understand this new enhancement 
structure observed in the baryonic B decay. In this work, we 
propose a novel approach to explain the BaBar's observation 
of X e (3250)°, which can be naturally explained as a molec- 
ular hadron composed of a charmed meson D ( *,(2400) and a 
nucleon N. In the following, we illustrate why the explana- 
tion of the D*,(2400)A^ molecular hadron is reasonable for the 
observed X c (3250)° in detail. 

A molecular explanation requires that the mass of the ob- 
served state should be below and close to the sum of masses of 
its components. The X f (3250) newly observed by BaBar just 
meets this necessary condition because the mass is near the 
threshold of D*>(2400) and N. For the D* Q (2400)N molecular 
system, the flavor wave function is written as 



|X C (3250)°> = |D ( *(2400)°«> 
i | X (3250)+) = l D S( 2400 > + »H D S( 240 °)>> 

|X e (3250) ++ ) = 



[Z}*,(2400?» 



(1) 



: |F e (3250) + ) 



|D ( *(2400) + n) + \D*(2400)°p) 



whose expressions correspond to isotriplet and isosinglet, re- 
spectively, where |X C (3250)°) is the flavor wave function of 
X c (3250)°, etc. Under the assignment of the D*,(2400)A^ 
molecular state, we further deduce the quantum number of 
X c (3250) as I(J P ) = 
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With the neutral X f (3250)° as an example, we first carry 
out the calculation of the binding energy for the D* Q {24QQ)N 
molecular system. According to the mass values of the 
D*(2400)° and neutron n listed in Particle Data Group (PDG) 

[2], we can obtain the binding energy of X c (3250)° as 13 

MeV, with Mjy ( 24oo)» =2318 MeV and M„ = 940 MeV. Then 



we have M D . (24 oo)° + M„ = 3258 MeV > M s 



(3250) 



In the fol- 



lowing, we will examine whether it is reasonable to explain 
X c (3250)° as the n£>*(2400)° molecular hadron. To deduce 
the effective potential between D* | (2400)° and «, we adopt the 
effective Lagrangian of the light mesons interacting with the 
charmed meson D*,(2400) or nucleon, i.e., 

£ mD >D> = -i^Dl b Di(2iv-Y ba ) + 2CcrDl )a Di, (3) 



inNN 



- yf2g VNN N b |y + ^d") Y^ ba N a + gcrNN NcrN 

(4) 

with v = (1,0), cr^ = Uy^y - /YO and the vector matrix V 



— PL 

V2 



P 



The index m in -Cmfljflj and -CmAW denotes the light meson. 
The coupling constants involved in this work are given by B' = 



I, 



0.76 [3, 4], g 2 pNN /4n = 0.84, gl NN /4n = 20, g v 



(2) m p lf n = 5.8, gi NN l4n = 5.69 and k = 6.1(0) forp(w) [5-8], 
respectively, where we follow the convention of the signs of 
coupling constants in Refs. [6-8]. The mass of the exchanged 
sigma is taken as 660 MeV. 

The effective potential of X c (3250)° in the coordinate space 
is given by 



'VTotaiM = i%(r) + i%,(r)+'V .(r) 



(5) 



with ^VV(r) = -2B'g V NNgvY(A,m v ,r) and ^(r) = 
^g'g.go-iw Y(A, mo-, r), where the Y(A, m, r) function is de- 
fined as 



1 

Y(A, m, r) = 

47rr 



-mr -Ar 

? - e 



A 1 



2A 



2 



By solving the Schrodinger equation with the effective po- 
tential obtained in Eq. (5), the dependence on the cutoff A 
of the bound state solution for X f (3250)° is shown in Fig. 1. 
When A = 1 .23 GeV, the theoretical result of the binding en- 
ergy for X c (3250)° is consistent with the experimental data 
[1]. Usually, in the one-boson exchange model the general 
criteria of forming a molecular state is that we can obtain neg- 
ative binding energy of this system and the corresponding cut- 
off A should be around 1 GeV. Thus, by our calculation, we 
can conclude that D*(2400) and nucleon can form a loosely 
bound state with the small binding energy since the adopted 
A value is close to 1 GeV. 
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FIG. 1: (color online). The obtained bound state solution (binding 
energy E and root-mean-square r RMS ) of X c (3250)° dependent on 
cutoff A. Here, the blue dashed line corresponds to £ = -13 MeV 
derived by the mass difference M Xl - M D « - M„. In addition, we also 
present the variation of the total effective potential and the subpoten- 
tials in Eq. (5) in r. 

As shown in Fig. 1, the potential of the Do(2400)°« molec- 
ular state is mainly from p exchange, where the n exchange 
is absent. It is well known that the vector meson exchanges 
provide short-range interaction. Usually, the necessary condi- 
tion of forming a bound state requires that the mass of the ex- 
changed meson is larger than the widths of the components of 
this bound state. We notice that the p meson mass (m p = 770 
MeV) is larger than the width of D (2400)° (r ~ 267 MeV), 
which basically satisfies the necessary condition of forming 
a molecular state. Considering the above reasons, Do(2400)° 
and n can interact with each other before Do(2400)° decays 
into other final states. This is why we still expect to discuss 
the molecular configuration to X c (3250)° although D (2400)° 
is a very broad meson. 

Since the total width of a state is mainly determined by 
its dominant decay, studying the dominant decay mode of 
X c (3250)° provides the important information of its total de- 
cay width. What is more important is that this study can be as 
a critical test of the nD ( *,(2400) molecular state assignment to 
X c (3250)°. 

As the «D*(2400) molecular state, X f (3250)° decays into 



£>;(2400)°. Then, decays Z c (3250)° -> D + '°n-°n occur 
via the intermediate Dj* ) (2400)° as shown in Eq.(6). Since 
the branching ratio of D*(2400)° — > Dn is almost 100%, 
X r (3250)° — > D +X) n~ f) n is the dominant mode. For this pro- 
cess, the differential decay width is written as 




{2n?2e x (27!-) 3 2e 2 (2n) 3 2e 3 



(6) 



with / = 1/2, where g,(e,) (i = 1,2, 3) denotes the momentum 
(energy) of final states. The decay amplitude M is expressed 
as 



M = 



^(^(3250)° -> D* Q n)Jl(D* -> Dn) 
q 2 -Mi, +iM D .Y D . 

1 D" u a u a 



(7) 



where M Xc(32 50) = -»(3250)° -» D* n) and ^ = f^D* -» 
Dn) describe the interactions X c (3250)° «D*(2400)° and 
D;(2400) -h. D + '°n-'°, respectively. M D; and T D . g denote 
the mass and width of the charmed meson D()(2400)°, re- 
spectively, q is the four momentum carried by the off-shell 
D;(2400)°. 

By the convariant spectator theory (CST), we can describe 
the collapse of X c (3250)° into the on-shell n and the D*(2400), 
where the vertex function |F) satisfies the relation 



W) = vc\r), 



(8) 



which is obtained by the Gross equation. In Fig. 2, we 
present the diagrammatic representation of the Gross equa- 
tion for the vertex function F. We further obtain the Gross 



P/2 + p 



P/2 + k 




P/2-p 



P/2 - k 



FIG. 2: (Color online.) Diagrammatic representation of the Gross 
equation for the vertex function T. Here, n and £>q(2400)° are marked 
by indexes 1 and 2, respectively, x denotes n being on-shell. 



equation, which is dependent on three-momentum p as 

d 3 k 



r(p) 



r d'k 
"J On) 



■<V(p, k, W)G(k, W) T(k) 



(9) 



D + '°n- 



Here, X c (3250) first falls apart into n and 



with F(p) = Y(p)ur(W). We use the convention liu = 2m 
here and hereafter, where we perform the integration over k°. 
P = (W, 0) denotes the four-momentum of the «D*(2400)° 
system, p - (pi - pi) 12 and k — {k\ - ^)/2 are the relative 
momenta as depicted in Fig. 2. Just because neutron is on the 



3 



mass shell, one gets k = (ko, fe), p = (po, p), Icq = £i(&) — 5 W, 

= j + k 2 . 'Vip, fe, W) is the interaction kernel with 
neutron on the mass shell. The two-body Green's function 
G(fc, W) in Eq. (9) is expressed as 



G(fc, W) 



-J—u®(k)u®(k) —L- ( l - 

2£i(fc) 2E 2 (k)\E 2 (k)- E x (k) + W 



1 



E 2 (k) + E l (k)- Wj 
The normalization of the vertex T(p) requires [9, 10], 



(10) 



—E-^(p)—[G(p,W)]T(p). 

(2nY 



(11) 

The wave functions for the nDg(2400) u bound state are 

1 r 



«a;(p) = 



1 



V(2tt) 3 2W \ 2E l (p)2E 2 (p) E 2 (p) + E 1 (p)-W 

(12) 

1 u^{pt(p)u R {W) 



V(2tt) 3 2W \ 2E l (p)2E 2 (p) E 2 (p) -Ei(p) + W 

(13) 

With the above preparation, the integral equations are writ- 
ten as 

[E 2 (p) + E l (p)-W] #(p) 
cfik 



[V r Ap, fe, W)^(fc) + V rr> (p, fe, W)^(fe)], 



/ (2tt)3 
[E 2 (p) - E!(p) + w] <a;(p) 

<f 3 A: 



-/ 



(2^ 



[v^(p, fe, wo^(fe) + y„-(p, k, w)<e;,(k)i 



where the expression of potential V rr >{p, fe, W) is 
V,Ap,k,W) = 



1 



1 



1 



1 



V2£i(p) \ 2E 2 {p) \ 2E x (k) \2E 2 (k) 

xu (r \p) r V(p,k,W)M (lJ \k). (14) 

After taking the nonrelativized approximation [11-13] and 
the Fourier transformation, we get the integral equations in the 
coordinate space 



+ ej (/,+(>) = -[V(r) + V(r)F(r)V(r)]ifr + (r)(l5) 

rir) = -[F(r)V(rM + (r), (16) 

where F(r) = [m2 - m\ + W + V(r)]~ l and E\{p) + E 2 (p) - 
W ~ —e + — with the reduced mass fi and the binding energy 
e = W — mi — m 2 . 

For the loosely bound system discussed in this work, it is 
reasonable to assume mi t 2 » (V). Then, we have 



-^ + eJ(A + (r) = -V(r)iA + (r), 



ifT(r) = 0, 



(17) 
(18) 



where Eq. (17) corresponds to the Schrodinger equation. We 
need to specify that V(r) is the total effective potential of the 
«Z)q(2400) molecular state shown in Eq. (5). Thus, by solving 
the Schrodinger equation with the deduced effective potential 
for the nD( } (2400) molecular state, we obtain if/ + (r). By the 
Fourier transformation, we get the wave function 4i + (p) in the 
momentum space 



<A + (P) 



(27T) 3 / 2 /' 



d 3 re- ipr ^ + (r), 



(19) 



which satisfies the normalization condition J cPp |<A + (p)| 2 = 1 
required by Eq. (11), where p denotes the relative momentum. 

With the calculated if/ + (p) and Eq. (12), we obtain the ver- 
tex F(p), which directly corresponds to JHx c (3250) in Eq. (7). 
By the CST, the calculated wave function is related to the ver- 
tex of the X c (3250)° collapse into «D;(2400)°. Thus, we find 
the relation 

-^X r (3250) 



V2W(2tt) 3 ^2E n {q) 2£ D j (240 0)(g) <A, + (<7) 
-W + E„(q) - E D > o(2 4oo)(q) 
In addition, the amplitude Jio' reads as 



(20) 



DJ°(2400) 



= Hg^M^oMo, (21) 



where the isospin factor / is taken as 1 and 1/ y2 for decays 
D^(2400) -> D + n- and D*(2400) -> D°n°, respectively. The 
coupling constant g„ is determined by the decay width {Vd- = 
267 + 40 MeV [2]), i.e., 



r(D ( *°(2400)° 



D + -V'°) 



m d Ap\ 

8ttM d . 



(22) 



with p being the three-momentum of the daughter meson in 
the rest frame of the D*°(2400) meson, where the Dn channel 
contributes most to the total width of 0^(2400) [14, 15]. 

In Table I, for several typical values of A we give the de- 
cay width of X f (3250)° — > nDn, which is the dominant de- 
cay mode of Z c (3250)° but does not strongly dependent on 
A. We compare the theoretical value of the decay width 
X c (3250)° -> nDn with the BaBar's data [1], which shows 
that the total width of X e (3250)° under the assignment of the 
nD ( *j(2400)° molecular state is comparable with the BaBar's 
measurement [1]. The study of the dominant decay channels 
of X c (3250)° also supports the «D*(2400)° molecular state ex- 
planation for the observed X c (3250)°. 

Apart from calculating the dominant decay width of 
X c (3250)°, we also calculate the corresponding line shape of 
the pion spectrum of the X c (3250)° — * nD + n~ process with 
the help of the CERNLIB program FOWL (see Fig. 3). 

Since X f (3250)° was observed in the Y.^ + 7i~n~ invariant 
mass spectrum, in the following analysis we illustrate how 
this process can happen under the «D*(2400)° molecular state 
assignment. Due to this molecular picture, X e (3250)° first dis- 
associate into off-shell n and D ( *j(2400) , which then transits 



4 



TABLE I: The decay width of X c (3250)° -> nDn dependent on A 
and the comparison of our result with the BaBar result. Here, we also 
calculate the binding energy E and the decay width T of 5 / c (3250) + , 
which is the isoscalar Dg(240Q)N molecular state defined in Eq. (2) 
as the partner of X c (3250). A, E and T are in units of GeV, MeV 
and MeV, respectively. The effective potential of y c (3250) + can be 
easily obtained by replacing the factor in front of "Vpir) in Eq. (5), 
i.e., 1/2 -> -3/2. 
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FIG. 3: (color online). The Dalitz plot analysis and the pion spectrum 
for the X c (3250)° -> n£> + 7T decay. 



into the final states E+ + 7r~7r~ by exchanging the D + meson 
(see Fig. 4 (a)-(b) for more details). Thus, we can naturally 
explain why X c (3250)° is reported in the l^. + n~n~ invariant 
mass spectrum. 




(a) (b) (c) (d) 




(e) (f) (g) (h) 



FIG. 4: (color online). The X r (3250)° -> decay and other 

possible decays of X ( .(3250)°. 

We also qualitatively discuss other possible decay modes of 
X c (3250)°, which are shown in Fig. 4 with the corresponding 
schematic diagrams. Here, the X c (3250)° decays into D°n, 



D*°n, A+7r~, A+p - , 2p7r~, £ ( tp~ occur via these triangle dia- 
grams listed in Fig. 4 (c)-(h). We notice that these neutrons in 
the final states of the X c (3250)° -> D°n, D*°n decays are neu- 
tral. Thus experimental search for these two decay modes of 
X c (3250)° is difficult. Besides decaying into £> w °n, X e (3250)° 
can decay into a charmed baryon plus a light meson such that 
pion and p. Searching for X f (3250)° by these predicted decay 
channels will be an interesting research topic. Due to the sub- 
sequent decays E e — » A e + n and p — » 27r, finally the processes 
of X e (3250)° into A+p~, £+7T, E+p~ are related to these final 
states of a A c baryon plus multipion. Comparing with these 
processes, X c (3250)° — > A.*n~ is a typical two-body decay. 
We suggest experiment to carry out the search for X e (3250)° 
by its A*?T" invariant mass spectrum. 

If X c (3250)° is a D (2400)°n molecular state, its dominant 
decay mode is Dnn. In addition, we also listed its other de- 
cay modes (see Fig. 4), which occur via hadronic loop. Thus, 
these decay modes are suppressed compared with its Dnn de- 
cay. Although Dnn is its dominant decay mode, the final states 
of X(3250) — > Dnn contain neutral neutron, which makes that 
it is difficult to find its Dnn decay mode in experiment. It is the 
reason why experimental firstly reported X(3250) in £+ + 7r ~n~ . 

In this work we also study the possible existence of the 
isoscalar partner of X f (3250), which is named as F c (3250) + 
whose flavor wave function is given in Eq. (2). We find that 
the obtained A value corresponding to its bound state solution 
is around 3.3 GeV, which is larger than 1 GeV. Thus, accord- 
ing to the criteria (A ~ 1 GeV), we can conclude that it is 
impossible to form F t ,(3250) + molecular state. 

Now let us draw a brief conclusion. Being stimulated by the 
recent observation of an enhancement structure X c (3250)° in 
the £^ + 7T~7r~ invariant mass spectrum of B — > Y^ + pn~n~ [1], 
we find that X e (3250)° can be well explained as the isovec- 
tor «Dq(2400)° molecular hadron, which is supported by both 
the analysis of the mass spectrum and the study of its domi- 
nant decay channel. Furthermore, the observed X c (3250)° — » 
+ tCtC can be described reasonably under this picture. In 
addition, we also mention several other possible decay modes 
of X c (3250)°, which can be studied in future experiments. We 
expect the contributions from BaBar, Belle, LHCb, and forth- 
coming Bellell, SuperB, which are the ideal places to further 
investigate the observed X e (3250)° in the B decay. 
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